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In this article two-loop QED corrections to the muon decay and corrections of order 
a 2 to the semileptonic decay of the bottom quark are considered. We compute the 
imaginary part of the four-loop diagrams contributing to the corresponding fermion 
D \ propagator in the limit of small external momentum. The on-shell condition is ob- 

tained with the help of a conformal mapping and Pade approximation. Via this 
method we confirm the existing results by an independent calculation. 

X 

1 Introduction 

The Fermi coupling constant, Gp, constitutes together with the electromagnetic coupling 
constant and the mass of the Z boson the most precise input parameters of the Standard 
Model of elementary particle physics. Gp is denned through the muon lifetime, and the 
decay of the muon, as a purely leptonic process, is rather clean — both experimentally 
and theoretically. The one-loop corrections of order a were computed more than 40 years 
ago H] whereas only recently the two- loop corrections of order a 2 have been evaluated || . 
The large gap in time shows that this calculation is highly non-trivial. The inclusion of the 
two-loop terms removed the relative theoretical error of 1.5 x which was an estimate 
on the size of the missing corrections. The remaining error on Gp now reads 0.9 x 10~ 5 and 
is of pure experimental nature. Upcoming experiments will further improve the accuracy of 
the muon lifetime measurement and therefore the O (a 2 ) corrections to the muon decay are 
very important and constitute a crucial ingredient from the theoretical side. These facts 
make it desirable to have an independent check on the correctness of the 0(a 2 ) result. We 
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also want to mention that in a recent article || optimization methods have been used in 
order to estimate the coefficient of order a 3 . 

In view of the upcoming B physics experiments the evaluation of quantum corrections to 
B meson properties have become topical. In particular it is possible to use the semileptonic 
decay rate of the bottom quark in order to determine the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix elements with quite some accuracy. In Q an approximate expression for 
the 0(a 2 ) corrections of T(b — > cW{) has been obtained where a non- vanishing charm quark 
mass has been included. In the decay b — > uev e the mass of the u quark can be neglected 
which reduces the calculation to the situation given in the muon decay. The only additional 
diagrams are those which arise due to the non-abelian structure of QCD. The results of 
order a 2 have been obtained in |J. The corrections proportional to the number of light 
quarks have already been computed in [|]. 

In this letter we confirm through an independent calculation the results of order a 2 to 
the muon decay || and the one of order a 2 to the semileptonic bottom quark decay . In 
the next section our method is discussed, in Section | the results are presented. 



2 Method and Notation 

In [H, U the imaginary part of the muon propagator is computed up to the four-loop level. 
Recurrence relations based on the integration-by-parts technique |7[] are used in order to 
reduce the integrals to be evaluated to a minimal set — so-called master integrals. They 
are finally evaluated by computing expansions in the ratio of external momentum and 
internal muon mass. It is possible to take the on-shell limit and perform the infinite sum 
which finally leads to an exact result for the integrals. For concise reviews of expansion 
methods see e.g. || ||. 

In contrast to P, [5| our approach is based on an expansion of the full fermion propagator 
in the limit 

M 2 > q 2 , (1) 

where q is the external momentum and M is the propagator mass of the muon and bot- 
tom quark, respectively^. Throughout the whole paper we will neglect effects induced by 
the non-vanishing electron and up quark mass. The on-shell limit q 2 — > M 2 will be per- 
formed afterwards with the help of Pade approximations. This, of course, only provides an 
approximation to the exact result. However, the integrals to be evaluated are simplified 
considerably. We will demonstrate that the accuracy obtained with our method is sufficient 
to check the existing results and enables the same reduction of the theoretical error on Gp- 
The notation is essentially adopted from flQfl , where corrections of 0(a 2 ) to the decay 



t — > Wb have been computed. For completeness we briefly repeat in this paper the main 
formulae. The decay rate — both for the muon and the bottom quark — can be written 



lr There are, of course, diagrams that do not contain an internal propagator of mass M. These diagrams 
are computed without any expansion. 
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M is the on-shell mass. TP S and Tj V represent the scalar and vector part of the corresponding 
fermion propagator, respectively. They are functions of the external momentum q and the 
bare mass mP of the fermion under consideration. In our case they further depend on the 
bare electromagnetic coupling a and the strong coupling constant a°, respectively, and 
are proportional to the square of the Fermi coupling constant, G 2 F . 
The mass renormalization constant, Zr~, entering in 



can be extracted from [11 



In contrast, Z° s has to be evaluated in the limit M 2 3> q 2 since the handling of Z® is 
determined by the computation of the fermion propagator. As we are only interested in 
the imaginary part and furthermore consider only QED or QCD corrections to the leading 
order term the quantity Z® s has to be known up to two loops only. The result can be 
taken from [ nj . The renormalization of a and a s proceeds in the usual way where we have 
chosen to renormalize also the electromagnetic coupling in a first step in the MS scheme. 

In order to get reliable results it is necessary to compute as many terms as possible 
in the expansion parameter z. Subsequently a Pade approximation is applied which is 
at length described in flOfl . We just want to mention that before the Pade procedure a 



conformal mapping can be used which maps the complex z-plane into the interiour of the 
unit circle. Following Ref. |10| we denote those results by tu-Pades and the ones obtained 
without conformal mapping by z-Pades. 

Some Pade approximants develop poles inside the unit circle (\z\ < 1 and \u\ < 1, 
respectively) in conflict with the analycity of the exact result. In general we will discard 
such numbers in the following. In some cases, however, the pole coincides with a zero of the 
numerator up to several digits accuracy, and these Pade approximations will be included 
in our sample. To be precise: in addition to the Pade results without any poles inside 
the unit circle, we will use the ones where the poles are accompanied by zeros within a 
circle of radius 0.01, and the distance between the pole and the physically relevant point 
q 2 /M 2 = 1 is larger than 0.1. 

The central values and the estimated uncertainty will be extracted from Pade results 
[m/n] with m+n not too small and \m— n\ < 2. The central value is obtained by averaging 
the Pade results and the uncertainty is given by the maximal deviation from the central 
value. 

It is convenient to parameterize the radiative corrections for the semileptonic bottom 
quark decay in the following form: 
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Figure 1: Sample diagrams for the fermion self energy. In the case of the bottom quark 
(thick line) decay two of the thin lines represent the lepton pair and the third one the up 
quark. All (one-particle-irreducible) diagrams involving the coupling of the gluons (loopy 
lines) to the up and bottom quark have to be taken into account. The modifications in the 
case of the muon decay are obvious. 
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(2) = C 2 F A^\ + C A C F A^ A + C F T ni A^+C F TA^ F , 



(5) 



n 1 



G|M b 5 |Kfe| 2 /(1927r 3 ). For QCD the colour factors are given by C F = 4/3, 



with 

Ca = 3, and T = 1/2. ni is the number of massless quark flavours and will be set to n\ = 4 
at the end. A\' A corresponds to the abelian part already present in QED, A\' NA represents 

(2) (2) 

the non-abelian contribution, and A b / and A b F denote the corrections involving a second 
fermion loop with massless and heavy quarks, respectively. In principle there is also a 
contribution involving a virtual top quark loop. It is, however, suppressed by M^/M t 2 and 
will thus be neglected. In Fig. [1] a representative diagram for each one of these four colour 
structures is pictured. In Eq. @ a s = a s (/i) is defined with five active flavours. The 
analytic result for A^ can be found in ||. 

An analogous expression to (Kf) can also be defined in the case of the muon decay 
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a(/i) represents the electromagnetic coupling in the MS 



scheme. Aff^ represents the purely photonic corrections whereas A^ e and Aff contain 
an additional electron and muon loop, respectively. The contribution involving a virtual r 
loop is not listed in (f|) as it is suppressed by M^/M 2 and almost four orders of magnitudes 
smaller than the other terms 



12|. 



Let us in the following describe our method used for the practical calculation in the 
case of the [i decay. The difference to the quark decay consists only in the transition from 
QED to QCD which increases the number of diagrams and makes it necessary to include 
the colour factors; the idea is, however, applicable in the very same way. 

Following common practice, we investigate the effective theory where the W boson is 
integrated out. The QED corrections to the resulting Fermi contact interaction were shown 
to be finite to all orders It is quite advantageous to perform a Fierz transformation 
which for a pure V — A theory has the consequence that afterwards the two neutrino 
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lines appear in the same fermion trace. Thus the QED corrections only affect the fermion 
trace involving the muon and the electron. This also provides some simplifications in 
the treatment to 75 since in the case of vanishing electron mass a fully anticommuting 
prescription can be used. 

As described above we consider the fermion two-point functions and compute the imag- 
inary part arising from the intermediate states with two neutrinos and the electron. As a 
consequence already for the Born result a two-loop diagram has to be considered. However, 
it turns out that the loop integration connected to the two neutrino lines can be performed 
immediately as it constitutes a massless two-point function. This is also the case after 
allowing for additional photonic corrections. As a result one encounters in the resulting 
diagram a propagator with one of the momenta raised to power e where D = 4 — 2e is 
the space-time dimension. This slightly increases the difficulty of the computation of the 
resulting diagrams. Especially for the order a 2 corrections, where the original four-loop 
diagrams are reduced to three-loop ones with non-integer powers of denominators, it is a 
priori not clear that these integrals can be solved analytically. However, it turns out that 
for the topologies needed in our case this is indeed possible. For the computation of the 
massless two-point functions we have used the package MINCER |TJ]]. Slight modifications 
enabled us to use this package also for the computation of the new type of integrals. 

The calculation is performed with the help of the package GEFICOM ]15[. It uses 



QGRAF ]nj for the generation of the diagrams and EXP [|17]] for the application of the 



asymptotic expansion procedures. For more technical details we refer to a recent review 



concerned with the automatic computation of Feynman diagrams |0 

The dispersive and absorbtive part of the fermion self energies are gauge dependent for 
q 2 7^ M 2 . Hence the dependence on the gauge parameter £ in Eq. (H) only drops out after 
summing infinitely many terms in the expansion around z = and setting z = 1. Since we 
are only dealing with a limited number of terms, our approximate results will still depend 
on the choice of £ even after taking z — *■ 1. It is clear that with extreme values of £ almost 
any number could be produced. Thus the question arises which value of £ should be taken 
in order to arrive at a reliable prediction for the decay rates. As the one-loop corrections 
can be evaluated for an arbitrary gauge parameter an extensive study can be performed 
and we can gain some hints for the choice of £ at order a 2 . 



3 Results 

Let us in a first step present the results for the muon decay and afterwards discuss the 
additional diagrams necessary for the QCD corrections to the bottom quark decay. 

The lowest order (Born) diagram can be computed directly. In this case the electron 
mass can be chosen different from zero and an expansion in M 2 /M 2 can be performed 
reproducing the exact result 

A® = 1-8— | - 12—^ In —| + 8— I 1 . 7) 

Ml M* Ml Ml Ml 1 ; 
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Table 1: z-Pade results for the 0{a) corrections for different choices of £. 



We want to demonstrate the power of our method for the order a corrections where 
four three-loop diagrams contribute. At this order we are able to evaluate a large number 
of moments which gives us a suggestion how many terms are necessary at 0(a 2 ) in order 
to obtain reliable results. Furthermore the computation can be performed for arbitrary 
gauge parameter which also provides some hints for the two-loop QED corrections. 

Applying the asymptotic expansion in the limit M 2 ^> q 2 to the four three-loop dia- 
grams contributing to the 0{a) correction leads to the following result for the first nine 
expansion terms 
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In Tab. |T] results for the 0{a) coefficient can be found where the Pade approximation is 
performed in the variable = q 2 /M 2 . Furthermore the gauge parameter defined through 
the photon propagator i(—g^ + ^q ll q 1J /q 2 )/(q 2 + ie) is varied between^ = —2 and £ = +2. 
Pade results which develop poles for |^| < 1 are in general represented by a dash. If an 
approximate cancellation with a zero from the numerator takes place (see the discussion 
above), they are marked by a star (*). The comparison with the exact result shows that 
for all values of £ reasonable agreement is found. However, there is a clear preference for 
£ = where the agreement with the exact result is best[]. Thus we will adopt this value 

2 Despite the fact that for £ > 1 the generating functional is in principle not defined we decided to 
choose this range for the gauge parameter. 

3 It seems that the ^-dependent terms of A$exp follow from the construction rule 25/48— Y!m=i z vl n l ( n + 
4) which for z il = \ indeed gives zero. 
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Table 2: Pade results for the 0(a) corrections (z- and cu-Pades). 

for the 0(a 2 ) calculation. 

In Tab. § the gauge parameter is fixed to £ = and in addition to the z-Pades also the 
cu-Pades are listed. With the inclusion of more moments the approximation to the exact 
result improves. Taking only those results into account where eight or nine input terms 
enter the following result for the order a correction can be deduced 

A® = -1.80(1). (9) 

Here the notation 1.80(1) = 1.80 ± 0.01 has been adopted. The excellent agreement with 
the exact result quoted in Tab. || encourages the use of our method. 

At order a 2 only the first eight moments are at hand. Using only seven and eight input 
terms changes the numbers of Eq. (||) to 

= -1.80(2), (10) 

which is still sufficiently accurate. 

Let us now move on to the two-loop QED corrections. Altogether 44 four-loop diagrams 
contribute. The application of the asymptotic expansion in the limit M 2 ^> q 2 leads to 72 
sub- and cosub-diagrams, which have to be evaluated. The analytical expressions obtained 
from the asymptotic expansion are quite lengthy. Thus we refrain from listing the results 
explicitly and present them only in numerical form. 

Our results at 0(a 2 ) are summarized in Tab. |3] where the scale fi 2 = M 2 has been 
adopted^. For comparison in the last line the numbers presented in Ref. [0] are listed. The 
results one obtains using the Pade approximants computed with seven and eight input 
terms read 

4 2 ) 7 = 3.5(4), 
45 = 3 -2(6), 

4 2 i = -0-0364(1). (11) 

4 Thc ln/i 2 /M2 terms can be reconstructed with the help of the (3 function governing the running of 
a(fi). 
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Table 3: Pade results for A^, and A<?\ 
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Table 4: Pade results for the corrections of 0(a 2 ) to the muon decay, A^\ 

It is remarkable that the central values agree well with the exact results which can be inter- 
preted as a sign that the presented error estimations are quite conservative. Furthermore 
it can be claimed that via our method we were able to confirm the results of 0. 

The error of A^' is particularly small, as the expansion in z converges very quickly. 
Note that in this case all a>-Pades develop poles inside the unit circle. The z-Pades are, 
however, very stable. A similar behaviour has been found for the analogous contribution 
to the decay of the top quark into a W boson and a bottom quark . 

A prediction for the decay rate of the muon up to order a 2 can be in principle obtained 
by summing the individual terms of ([TT|) . This would, however, significantly overestimate 
the error. It is more promising to add in a first step the moments of the single contributions 
and perform the Pade procedure for the sum. The corresponding results are shown in Tab. |] 
which finally lead to 

4 2) = 6 - 5 ( 7 )- (12) 

The deviation of the central value from the exact result of 6.743 is less than 3% and well 
covered by the extracted error of roughly 10%. Thus the sole knowledge of our results 
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Table 5: Pade results for the non-abelian part of 0(a 2 ), A[' NA . 

would also reduce the theoretical error on Gf as mentioned in the Introduction. Using the 
results presented in this paper the decay rate of the muon reads 

r(/x -> z/ M ez/ e ) = rj 

As already noted in H] the numerical coefficient in front of the second order corrections 
becomes very small if one uses the on-shell scheme for the definition of the coupling constant 
a. Then the MS coupling is given by a(M M ) = a(l + a/37rln(M^/M ( 2 )) and there is 
an accidental cancellation between the constant and the logarithm in the second order 
corrections. 

Let us now turn to the semileptonic decay of the bottom quark. The Born and one- 
loop corrections can, of course, be taken from the muon decay. In particular we have 
= A)®' and = A^' . As far as the two- loop terms are concerned only the non- 

(2) 

abelian contribution, A K h NA) has to be computed in addition. The other colour structures 
are related to the expressions occurring in the muon decay rate through Aj/ A = A@), 
AfJ = A^ e , and A^ F = A^L with obvious replacements of the masses. 
In Tab. |5] the results for A^ NA can be found. We infer 

A® NA = -8.8(4), (14) 

where /i 2 = M 2 has been chosen. The central value is again in very good agreement with 
the exact result and agrees within the error estimate of 5%. 

In order to get predictions for A b we again add the moments in a first step and perform 
the Pade approximations afterwards. From the results listed in Tab. || we deduce 

4 2) = -21.1(6). (15) 

This number is in good agreement with the one stated in |J. The error is quite small 
and amounts only to roughly 3%. The semileptonic decay rate of the bottom quark finally 



0.9998 - 1.810 



71 



+ 6.5(7) 




+ . . . 



(13) 



9 



in nnt, 

J.1J. kj LI u 


P.A. 


z 


CO 


6 


[3/2] 


-20.587 


-21.835 


6 


[2/3] 


-20.592 


-22.252 


7 


[4/2] 


-20.744 


-21.159 


7 


[3/3] 


-20.836 




7 


[2/4] 


-20.757 


-21.649 


8 


[4/3] 


-20.964 


-21.213 


8 


[3/4] 


-20.965 


-21.608 


exact: 


-21.296 



Table 6: Pade results for the corrections of 0(a 2 s ) to b — > uev e , Af\ 
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4 Conclusions 

In this paper the two-loop QED corrections to the decay rate of the muon have been evalu- 
ated. A new method has been used in order to confirm via an independent calculation the 
result of Ref. 0. From the muon decay the Fermi coupling constant, Gf, is determined 
which constitutes one of the basic input parameters of the Standard Model. Thus it is 
very important to have independent checks on such highly non-trivial computations. The 
inclusion of the new correction terms removes the theoretical error on Gf- After the addi- 
tional computation of the non-abelian diagrams the QCD corrections to the semileptonic 
decay rate of the bottom quark T(b — > ueu e ) are obtained. Again agreement with the 
literature [||] is found. 
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